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Abstract 

The problem of pseudo-gap formation in an electronic system, induced by the fluctuations of 
the order parameter is revisited. We make the observation that a large class of current theories 
are theoretically equivalent to averaging the Free energy of the pseudo-gap system over quenched- 
disordered distribution of the order parameter. We examine the cases of both infinite and finite 
correlation length, showing how the interplay of pseudo-gap formation and superconductivity can 
be treated in this approach. 
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I. INTRODUCTION 



The discovery of a marked suppression in the density of states of the underdoped cuprate 
superconductors above the critical temperature has led to a revival of interest in the idea of a 
pseudo-gap. This concept was originally introduced by Lee, Rice and Anderson |j] to explain 
the suppression of electron density of states associated with order parameter fluctuations 
near a charge density wave instability. According to one school of thought, a similar mech- 
anism may drive the formation of the pseudo-gap in the underdoped cuprates(see reviews 
3, El Many candidate fluctuating order parameters, such as antiferromagnetismjsj], 
charge density wave[(| and pre-formed Cooper pairs have been proposed 

Pseudo-gaps are likely to be a widespread feature of correlated electron systems lying 
close to an instability that gaps part of the Fermi surface, and they have been observed 
in a wide variety of nested and low dimensional strongly correlated electron materials, in- 
cluding one dimensional charge density wave systems]^], vanadium doped Chromium jjj, 
colossal magneto- resistance compounds [ic| and strontium-calcium and strontium-barium 
ruthenates|lll|. 

The fluctuation gap model (FGM), based on the idea of the pseudo-gap emergence due to 
fluctuating order parameter was pioneered in jlj]. Later on it was shown by Sadovskii [12I IlS^ 
how the problem of nested electrons moving in a critically fluctuating order parameter field 
could be solved exactly. Interest in this model was revived with the discovery of the pseudo- 
gap in cuprate, the physical origin of which „ S ti.l controversial Schma.ian et a. Q 
extended the ID model to higher dimensions and nonzero spin. Tchernyshyov [15|| later 
questioned the exact solvability of the model for finite correlation lengths and pointed out 
a previously unnoticed mistake in the Sadovskii result. Couple of years ago Millis and 
MonienQ] performed a careful study of different approximations, explaining for pseudo- 



ly oi 





gaps, also by Bartosch and Kopietz |17[ the exact numerical calculation of the density of 
states in FGM was carried out. They demonstrated that Sadovskii approach is in fact 

a good approbation and the term, omitted in his mode, produce negligible corrections. 
In a related development, Kopietz et al [19( presented the exactly solvable model of extended 
FGM, where fluctuations of both phase and amplitude of the order parameter have been 
taken into account. The important issues of crossover from Gaussian to non-Gaussian order 
parameter fluctuations in one-dimensional Peierls systems have been recently discussed by 
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A key idea of the Lee, Rice and Anderson paper pjis that order parameter fluctuations can 
be regarded as classical degrees of freedom with time-independent correlation functions. This 
basic assumption threads through most of the subsequent developments during past twenty 
years. This paper revisits this basic idea, making the observation that a large class of current 
theories are theoretically equivalent to averaging the Free energy of the pseudo-gap system 
over a quenched- disordered distribution of the order parameter. The intuitive equivalence 
between the slow critical fluctuations responsible for pseudo-gap formation and a quenched 
random potential was certainly known to several authors in the pastQ, 21, 13, 1^, I^J. In 
this paper, we vocalize this equivalence and set it in a formal framework. Using the idea 
we show how the Sadovskii model of pseudo-gap formation can be formulated in terms of 
a quenched average free energy. To illustrate the utility of this method, we develop a toy 
model for the interplay of pseudo-gap formation and weak impurity scattering with d-wave 
superconductivity. Our results may be of interest to recent experiments on the cuprate 
superconductors. 



II. CLASSICAL INTERACTIONS IMPLY QUENCHED DISORDER 

One of the central assumptions of a large class of pseudo-gap models, is that the order- 
parameter fluctuations which scatter electrons to form the pseudo-gap are so slow that they 
can be considered to be infinitely retarded. Such a state of affairs could come about in 
a variety of ways. The most pragmatic point-of-view, is that the characteristic relaxation 
timescale of the fluctuations top is not actually infinite, but may be treated as such because 
it is much greater than the inelastic scattering rate of the electrons r e 

r p » T e . 

The result of this condition, is that electrons perceive the order parameter as a frozen degree 
of freedom, that may be treated as a classical variable. 

We now need to consider the interactions between the electrons induced by these classical 
order parameter fluctuations. In a Feynman diagram, the net interaction induced by these 
fluctuations has a number of important features (Fig. 1.): 

• It transfers no energy, and is thus infinitely retarded. (Fig. 1 (a)) 
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• Polarization bubble insertions into the effective interaction are absent, since these have 
already been taken into account in forming the effective interaction (Fig. 1(b)). 

• Interactions between the classical fluctuations may be ignored. This is a useful, as- 
sumption, for it leads to non-interacting Gaussian fluctuations and permits the sum- 
mation of a large class of diagrams. Unfortunately, it can not really be justified near 
a critical point, where non-linear interactions between classical fluctuations are al- 
most always relevant. This is without doubt a major weakness of the current class of 
theories, and one we shall return to in our final discussion. 



(a) -W 2 S(q) 




FIG. 1: (a) Classical fluctuations induce an interaction V e ff = — W 2 S(q) on a single line, (b) 
Fermion loop insertions into the interaction line are neglected which means that the average over 
the classical fluctuations is a quenched average, (c) By carrying out a quenched average, we also 
eliminate interactions between fluctuations of the form shown here. 

The first feature can be incorporated into a model by introducing a coupling into the Hamil- 
tonian 

H[<j>] = H + Wj2p?<f>-?- (!) 

In this paper we will imagine that H is a Hamiltonian in which effects of interaction have 
been taken into account by a mean-field approximation- so that for instance, H might be 
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a BCS Hamiltonian with explicit pairing terms. In ([!} — X]fc c j: - c k ls ^ ne density of 
electrons at momentum q and ip^ is a purely classical field representing the order parameter 
fluctuations. This quantity is to be averaged over the distribution 

12 1 



P[4>] = exp 



2S(q)_ 

and it thus gives rise to an effective, and infinitely retarded (attractive) interaction of the 
form 

Veff(q) = -W 2 S(q) 

At first sight, it might appear sufficient average the partition function over P[<f>], i.e 

z = J P[<t>}z[<f>) 

Z[<j)] = Tr [ e ~ pm ] . (2) 

This is an annealed average. However, such a scheme allows fermion loop renormalizations 
and non-linear interactions to develop within the interaction lines as illustrated in Fig. 1 
(b). These terms are to be dropped, under the assumptions mentioned above. To ensure 
this, we must carry out a quenched average over the classical field, which means that we 
average the Free energy, rather than the partition function 

F = -T I P[0]logZ[0] (3) 



This procedure is implicit to the approach of Sadovskii, and subsequent developments, by 
the selection of Feynman diagrams containing a single Fermion line. One way to implement 
the quenched average is to use the iV — > trick, 

F = -Tlim| P[0]^[0f-1). (4) 

The term containing Z[(f)] can be considered as the partition function of N identical replicas 
in the frozen classical field <fi. Formally, this is accomplished by introducing iV replicas of 
the fermion fields, labeled by the index AG [1, N]. The partition function is then given by 

Z N [</>]=Ti[e-f )HW M] 

where 

H^]= £ [H^ + wJ2p^-X (5) 

X=1,JV \ a I 



is the Hamiltonian of N identical replicas coupled to a single classical field. Diagrammatic 
contributions to ^[0]^ containing n fermion loops, scale as 0(N n ), so that in the iV — > 
limit, the only terms surviving in the diagrammatic expansion of the Free energy, are those 
with a single fermion loop, as follows 

/ O(N) W \ 



AF 



-T lim 



1 

N 





V 




(6) 



where AF = Fj — Fq is the change in the Free energy due to turning on the interactions, 
Thus the only terms which survive the N — > limit are those with a single fermion line 
or 

loopQ- 

In this way, all RPA renormalizations and non-linear interactions between the 
fluctuation fields are eliminated, as follows, 



AF = — T 







(7) 



where the replica index has been eliminated from the fermion lines. The advantage of this 
approach is that we can now start to formulate a theory of the pseudo-gap in the language 
of an effective Free energy. This becomes particularly useful when considering the case of a 
superconductor with a pseudo-gap. 



III. EFFECTIVE ACTION IN THE PSEUDO-GAP MODEL WITH INFINITE 
CORRELATION LENGTH 

To illustrate the above approach, we consider a two-dimensional electronic system with 
a nested Fermi surface spanned by the commensurate wavevector Q = (tt, it) as illustrated 
in Fig. 2. Fluctuations of the underlying order parameter are governed by a correlation 
function of the form 

fluctuations com 



e the nested Fermi surfaces via the effective interaction V e ff = — W 2 S(q). 
It was shown |25| that in the limit £ —>■ 00 the model under consideration can be solved 
exactly and the effective interaction takes the oversimplified form 

Veff(q) = -W 2 {5(q x - n)5(q y - tt)}. (9) 



FIG. 2: Nested Fermi surface with commensurate fluctuations considered in the toy model of this 
paper. 

In this paper, we will compliment the earlier work[25] by showing how the diagrams for the 
Free energy may be completely summed in this limit. The class of models we shall discuss 
take the form 

H = H + Hj (10) 

where Hq is a Hamiltonian which does not scatter between the nested Fermi surfaces and 
Hi contains the coupling to the density fluctuation modes. We begin by considering the 
simple model where 

H o = J2 e k c L c k« 

k 

k.q 

where (p^ is a classical field described by quenched average over the probability distribution 
P[4>] given above. If we denote the effective action by a wavy line, then the Feynman 
diagrams for the expansion of the Free energy are 
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(12) 

Here we have implicitly taken the iV — > limit of (@J), thereby eliminating all polarization 
renormalizations of the interaction lines. We have used the pseudo-momentum labels " — " 
and " + " to denote electrons lying on the "left" or "right" -hand side of the nested Fermi 
surface. An essential feature of our derivation is that the pseudo-momentum alternates 
between scattering events. (For incommensurate fluctuations, this requirement eliminates 
some of the diagrams such as the fourth diagram shown above). In the limit of an infinite 
correlation length, each interaction line can effectively be replaced by a static modulated 
scattering potential which scatters at the nesting wavevectors. Consider the simplest W 2 
diagram. In the limit of infinite correlation length, the interaction can be replaced by a 
delta function in momentum, so that 




W 

(13) 

where we have represented the vertex of the static scattering potential by a "cross" . When 
we come to consider higher order diagrams, each term of a given order gives rise to the same 
contribution, however, we must be careful to take into account symmetry factors. Let the 
number of closed diagrams of order W 2n be a n , then 




Around the perimeter of each diagram of this order, there are n propagators of both "— " 
and "+" pseudo-momentum. If we differentiate the sum of all 2nth order diagrams with 
respect to G^, we generate the 2nth order self-energy diagram, as follows 
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(15) 



But in the limit of infinite correlation length, this is equal to 



SGj. 





\w 


k+Q A 1 




wi 2 


J k+Q 


k V 3 • 


2n-l / 









X— X- 



-X X 



(16) 



Now the number of 2n th order self energy diagrams is given by the number of ways of 
connecting n interaction propagators to 2n interaction vertices so that 




> = (2n- 1)! 



x— *- 



(17) 



enabling us to identify a n = (2n — 1)!!. Using the relation 

where we denote (x — |)! = r(x + |), we can absorb the factor (2n — 1)!! into an integral 
over a Gaussian distribution of gap sizes as follows 



1 f 00 

W^ 2n (2n- 1)!! = \ d(e-< 2 / 2 ((W) 2n , 



in other words 



(18) 




d(e 



-C 2 /2 
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(19) 

Remarkably- all that is left of the average over the classical fields is a residual distribution 
over pseudo-gap sizes. In other words, in the limit of an infinite correlation length, the 
quenched average over the distribution of scattering potentials is replaced by a single average 
over a Gaussian distribution of scattering strengths. The quenched averaged Free energy is 
then given by 




(20) 



where F a = — T^2k n l n [ — G ^n)] is the Free energy of the "non-interacting" system. 
Combining both terms, we then obtain 

1 



F = 



J 



d(e-^ 2 F[(, W] 



(21) 



where 



F[{,W] = -T J2 Trln 

n, A;g|bz 



-G-\k-Q/2,iu n ) (W 

(W -G-\k + Q/2,iu; ri 



(22) 



is the Free energy for a gap size (W. The most important feature of our derivation above, is 
the conservation of the pseudo-momentum associated with each nested Fermi surface. This 
derivation can be easily generalized to include pairing (which conserves momentum ) and a 
restricted class of disorder scattering which scatters between states on the same nested side 
of the Fermi surface. We now examine these two cases in detail 



A. Pairing in the presence of a pseudo-gap 

To consider pairing, let us now take 

IA 



2 



k k 

where 7^ = cos(k x ) — cos{k y ) is a d-wave order parameter whose nodes bisect the nested 
Fermi surfaces. We have included an additional term |A| 2 /g to take account of the BCS 
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decoupling of the pairing interaction. To make the model more interesting, we shall consider 
a pseudo-gap with d-wave symmetry, 

k,q 

where 7g = cos(fc x ) — cos(k y ) as before. Since the pairing does not change the momentum 
of particles, the reasoning for a normal state is quickly generalized by using the appropriate 
Nambu electron propagators. Let us first derive the free energy expression for this case. The 
mean -field Hamiltonian can be expressed in the following general form 

IAI 2 



kd\Bz 
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(24) 



where 



( C k+ Q/ 2 ] \ 



u -k-Q/2l 
V C -k+Q/2l 



(25) 



and 



^ e k+Q/2 ^k+Q/2 W k 



Ti 



\ 



(26) 



where W K = smK x — 
—W_ K , we can write 





— (ek+Q/2) —W-k 

W k Cfc-Q/2 Ak-Q/2 

-W-k ^k-Q/2 — e_(fc_Q/ 2 ) / 

sin«; y . Now since e K _ Q/2 = e_( K+Q / 2 ), A K _ Q/2 = -A K+Q/2 and W„ = 



V 



so that 





^ e K +Q/2 A k+ q/ 2 


w K ^ 




H = 


A k+ q/ 2 — e K+ Q/ 2 


+w K 






w K 


—£k+Q/2 —A k+ q/ 2 






V w« 


—A k+ q/2 e K+ g/ 2 J 




l {iuJ n )} 


= In det [H — iu n ] 


= H"l + 4+Q/2 + ^l+Q/2 + W 2 k 



(27) 



so that the the expression for the quenched-averaged Free energy takes the form 



V2n J 



d(e 



-C 2 /2 



-2T H^ + 4 +Q/ 2 + ^l +Q/ 2 + (CWk) 2 ] + — 



(28) 



11 



The gap-equation is provided by the condition that the derivative of the free energy is 
equal to zero, |£ = 0, i.e 

\ + -kr d ^ e A- 2T E ^tcttct^ ) = - (29) 

At the critical temperature T c , A = and we have the transition temperature equation 

9 V2i7-„ ^ „, te |,z aJ " + e * + C '^-0/2 J 

In this way, we have successfully eluded the diagrammatic derivation of the two-particle 
Green's function j2] and derived the main result in a compact and more general way. 



B. Impurity effect 

As a further application of these methods, let us now introduce disorder into the paired 
pseudo-gap model of the last section. The interplay of disorder and pseudo-gap formation 
on pairing is of particular interest in the context of underdoped cuprate superconductors. 
In a recent experimental paper, Tallon et al. have shown that the superconducting 
transition temperature of La — 241 and (Y, Ca) — 123 is much more rapidly suppressed by 
disorder in underdoped compounds, where a pseudo-gap is present. Tallon et al explained 
their results in terms of a phenomenological model involving unitary scattering taking place. 
We now show how the essential features of their discussion can also be obtained by treating 
the co-existence of pseudo-gap formation and weak non-magnetic impurity scattering. 

In the presence of impurity scattering we have to include the additional terms in our 
Hamiltonian of the form 

j,kk'cr 

In order to apply the methods of the last section, we make the important assumption that im- 
purity potential Vkk' does not scatter between different sides of the Fermi surface. Such scat- 
tering can still transfer electrons across the node, and thus remains severely pair-breaking. 
With this simplified assumption, the impurity scattering does not renormalize the strength 
of the pseudo-gap potential. 
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We can write now 

i poo 

F = ^= d(e~^ 2 Fi(,V kk ,}. (31) 

In this oversimplified model the effect of impurities reduces just to replacing of thin lines to 
there "thick" , impurity dressed equivalents in the graphic representation for the two-particle 
Green's function in our model. 

Using this result, we can generalize equation ([5U1) to obtain 

l = T T f f— ^ (32) 



g ^J-ooV^ Jo (2n)^l + el + CW^) 
where uo n is the usual renormalized Matsubara frequency 

u n = u n + Tsignuj n (33) 



T = 7^(0)1^, (34) 

and p is the impurity concentration. 

Notice that the momentum dependence of Wk-g/2 oc jk, so that in equation (JH2J), the 
functions 7(0) and W(<f>) share the same angular dependence around the Fermi surface. For 
simplicity we assume the d-wave symmetry of the pseudogap as well as of the superconduct- 
ing gap function 

W{<f>) = Wj{cf>), (35) 
where 7(0) = \/zcos(20) is the angular dependence of the gap around the Fermi surface 
EL 

After integrating over energy and carrying out the Matsubara sum, (Appendix A) we 
obtain the final equation for the critical temperature: 

^-r^>£^{'G)-G +E ^)}- (36 > 

where the normalized distribution function <$>(x) can be written in terms of the modified 
Bessel function K (x) as 

* w = 7h e '"" K ° (t) • 

The relative factor of " i" between the impurity and pseudogap term in this expression 
reflect the different mechanism by which impurity scattering and the pseudogap supress su- 
perconductivity. On the one hand, impurity scattering produces pair breaking, by scattering 
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particles across the nodes of the order parameter. By contrast, the pseudogap suppresses 
superconductivity by removing states from around the Fermi energy. 





FIG. 3: (a) Critical temperature of a superconductor, normalized by the critical temperature of 
a pure superconductor without pseudogap T c q, vs impurity scattering T in the presence of the 
pseudogap (W is an effective pseudogap width), (b) Critical temperature of a superconductor vs 
the effective width of the pseudogap for different values of impurity scattering rate r(both W and 
r are normalized by T c q ). These curves were all obtained from a numerical solution of equation 
(EH). 
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The theoretical results for the critical temperature dependence on the impurity scattering 
and the pseudogap, obtained by a numerical solution of Eq . (j3T?j) . can be presented in different 
ways. Fig. 3a, for example, shows the variation of the transition temperature with impurity 
scattering rate for a set of different values of the pseudogap. In the limit where the pseudogap 
vanishes, Eq. (J36|) naturally reverts to the standard Abrikosov-Gorkov mean- field theory. Fig. 
3b shows the dependence of the critical temperature on the pseudogap width for different 
ratios r = T/T cQ . These figures illustrate how both the pseudogap and impurity scattering 
suppress the superconductivity. 

To illustrate this approach we now use it to develop a phenomenology for the cuprate 
phase diagram (Fig. 4). We consider a model, in which the bare superconducting transition 
temperature is a linear function of doping x, vanishing at x — 0.3 and the pseudogap is a 
rapidly reducing function of x, which vanishes at x — 0.19. The effect of the pseudogap 
is to suppress superconductivity at low doping, leading to a maximum in the transition 
temperature in the vicinity of where the pseudogap goes to zero. When disorder is introduced 
the superconducting region is pushed to higher doping. 



IV. FINITE CORRELATION LENGTH 



We now return to the problem of pseudo-gap formation in the presence of order parameter 
fluctuations with finite correlation length. This problem was originally considered by Lee 
Rice and Anderson jl]. Following their approach, we consider Gaussian fluctuations of a 
static order parameter with finite correlation length and a real-space correlation function of 
the form 

< q*{x)V{x') >= \V\ 2 exp{-K\x - x'\)e i2pF{x - x '\ (37) 

where |\1/| 2 is the mean square fluctuation of the order parameter field and k is the inverse 
correlation length. 

Sadovskii jlj, [l| later extended the FGM model by considering the effect of multiple 
scattering off order parameter fluctuations. The effect of superconducting pairing was incor- 
porated into the model in recent papers [3(|. In the case of infinite correlation length, the 
Sadovskii yields an essentially exact solution in one dimension. The Sadovskii approach is 
only approximate at finite correlation length jlflj] . For example, Bartosch and Kopietz ^ 
have shown, that the Dyson singularity in the density of states, which exists for any finite 
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FIG. 4: Toy phase diagram of cuprates. Here W is the effective energy scale of the pseudogap, for 
simplicity we assume that this quantity follows a linear dependence on doping W = W max (l— x/x\) 
(x\ ~ 0.19), while we let the critical temperature of a pure superconductor T c q vanish at x ~ 0.3, 
since in most experiments superconductivity is not observed for doping value exceeding 0.4-^0.5. It 
is clearly seen now how the superconducting area shrinks while the scattering by normal impurities 
increases (T is impurity scattering rate (34), normalized by T c q). This phase diagram is surprisingly 



similar to the figure 5 of experimental paper 



29] 



value of £, is missed by Sadovskii's algorithm. Nevertheless, the Sadovskii approach seems 
to be justified in 2D case for certain topologies of the Fermi surface [^tJ. Nevertheless, the 
Sadovskii approach appears to provide a remarkably good approximation to the density of 



P,Q- 



states in a pseudogap system with finite correlation length 

To overcome some of the difficulties encountered in these early treatments, Kopietz has 
recently suggested a simpler alternative model, which exhibits the same static spatial cor- 
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relations originally considered by Lee Rice and Anderson, but which can be treated exactly. 
Kopietz's approach considers an order parameter of the form 



*(x) = Ae iqx (38) 

where the amplitude is a Gaussian random variable but the q— vector is determined by a 
Lorentzian distribution, as follows 

^^'°^ i?TFw (39) 

where we have chosen a one dimensional example, setting Q Q = 2pp. The order parameter 
correlation function is then 

(y*(x)V{x')} = J dAdAdqP(A, A, q)AAe iq{x - xl) = ^V«|z-zVQo(*-z') . (40) 

For the commensurate case this recovers the form (|H7jl first assumed by LRA. 

The Hamiltonian which describes the coupling between the order parameter fluctuations 
and the electrons is given by 

H(A, A,q)=J2 e kCkCk + J2^ AC k+q/2 Ck ~l/ 2 + Ac l- q /2 C k+q/2\- (41) 

k k 

Unlike the scattering described in (9), here we have a quenched distribution of critical 
scattering q— vectors, and it is this feature that gives rise to the finite correlation length. By 
resuming all diagrams associated with scattering off the fluctuating field, Kopietz was able 
to derive the spectral density for this model. 

We now show how the methods we have developed in the previous sections can be used 
to directly derive the Free energy for this model, from which secondary properties, such as 
the spectral Green function can also be derived. The starting point for our discussion is the 
Free energy in a given static configuration of the order parameter, given by 

F(A, A, q) = -T In Tr [ e -/»(^A<?)] ( 42 ) 

Following our earlier discussion, the Gaussian fluctuations over the static classical fields are 
now represented by the quenched averaged Free energy 

F = J dAdAdqP(A, A, q)F(A, A, q) (43) 
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The integrals over A, A and q can be carried out analytically (Appendix B) and yield 
the following result 



where u n = u n + (vFK/2)signu n . Thus the most important effect of the finite correlation 
length is to introduce an imaginary scattering term inside the Greens function. Notice that 
the distribution over gap sizes has become exponential (Rayleigh), rather than Gaussian, as 
it was in the case where Q Q = tt : this is a consequence of having assumed an incommensurate 
Q vector. Thus the essential physics of a finite correlation length manifests itself as an elastic 
impurity scattering potential. 

Various extensions of these arguments are possible, such as the extension to two dimen- 
sions, the use of a commensurate scattering potential and the introduction of pairing terms, 
but we shall not pursue them in this paper. 

V. DISCUSSION AND CONCLUSIONS 

The central idea of this paper, was to point out that a large class of treatments of the 
pseudogap problem may be formulated in terms of a quenched average over a distribution of 
static classical order parameter fluctuations. We have illustrated this method, showing how 
it can be used to construct the Free energy of the pseudo-gap system, from which mean-field 
equations can be derived. The case of infinite correlation length is easier to handle, but a 
simplified version of the case with finite correlation length has also been treated. 

The advantages of this method are that they permit us to begin with the Free energy 
functional, rather than working directly with Green functions. Thus we are able to introduce 
superconducting pairing into the formalism without having to reconsider a whole class of 
two-particle correlation functionsj^]. 

Despite the utility of the new method, there are a number of important questions that 
arise as to its complete validity. There are two obvious weaknesses in the current approach 

• The method ignores any frequency dependence in the fluctuations. For this reason, 
the method is only applicable at high temperatures. 

• A secondary weakness is the assumption that the classical modes are strictly non- 
interacting. With this widely used assumption, we have shown that the fluctuations 




(44) 
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behave as a source of quenched disorder. The equivalence between classical fluctuations 
and quenched disorder only holds when the modes are non-interacting and classical 

Both weaknesses are clearly issues of time scales. So long as the characteristic time-scale 
of the electrons is longer than the characteristic time scale of the fluctuations, then it is 
reasonable to treat them in the way we have outlined. If by contrast, we are interested 
in situations where the characteristic electron time scales become very large, such as the 
vicinity to a quantum critical point, then we might expect the methods used here to become 
invalid. 

Nevertheless, there are clearly a wide number of applications for this approach. One of 
the interesting possible applications are transport properties of a pseudogap system, and 
these are currently under active investigation. 
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VI. APPENDIX A 

The purpose of this appendix is to carry out the various integrations inside the gap 
equation 



(45) 



Q/2 



where we have denoted 



If the pairing is dominated by processes near the Fermi energy, we may replace the momen- 
tum sum by an energy integral, thus 
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so the gap equation becomes 

1 f°° r rift ( N max — 1 -, ^ 

- = 2nN(0)T c / d(P(() / ^ 7 (0) 2 { Y - \ . (47) 

where we have introduced an upper cutoff N max = D/2nT c into the Matsubara sum. To 
convert the Matsubara sum into a contour integral, we need to identify the function G(^)that 
has poles at z — iu n , (n G {0, N max — 1}). By using the identity 



W = -c + Efr —) 



we see that 



so that 



V" — ; = ^(z + N max ) -i/j(z), 



n=0 



°« = E — n = * (5 + -asr" J - *' (2 + (48) 

has the required properties. With this result, we can rewrite the Matsubara sum as 

Nmax 1 r rl? f \ 

2vr T c V - = 2ttT c / — - [iG(z)] 

— ; = (49) 

where the integral is taken anticlockwise around the positive imaginary axis. Notice that 
we have analytically extended iuj n — > z by replacing ^oj 2 + ( 2 Wg — > —iy/z 2 — ( 2 Wg. With 
this choice, there is a branch-cut running from z = — £W,9 to 2; = C^V Above and below 
the branch-cut, y/{x ± i5) 2 - ( 2 W 2 = ±i x /( 2 W 2 - x 2 . We may now distort the contour so 
that it runs clockwise around this branch-cut, to obtain 

N ™ a * 1 rCW e j I 

2nT Y - = / — - -.G[x]. (50) 

Along the branch-cut, we may replace 



A p + r + u\ f_D_\ 

^ \2 2nT J \2nT J 



so that inside the integral we can use 



/ D \ fl T + ix 

In 



^-nsfJ-^-sr 1 (51) 
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FIG. 5: Showing how contour around the poles of G(z) on the positive imaginary axis in equation 
(49) is distorted around the branch cut running between ±(Wg. 

At this point, the gap equation takes the form 

_l_ = I * (9) . r d(Pi0 1 r * 1 g[x\ ) . (52) 

We can further simplify this integral by reversing the order of innermost ( and x integrations. 
We do this as follows 

/ = r r —[...]= t/4 r° A* r \...\ m 

Aoo V 7 ^ J-CWe 7T V 7T J ~co J JM. 

Next we make the change of variable, w = ^ — ^fj^, so that the inner integral becomes 



I o /*oo foo — w 

,=^/- (54) 



Replacing the argument inside the integral by 

G[x] G[x] 



[■■■] 



we obtain 
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1 



V27T J J-oo 



dxe- x2/2 G(xW e ) / du 



l ^u(u + x 2 /2) 
The integral on the right hand side can be carried out analytically, 

'du . 6 ~ U =e t ' 2 K {l) 

y/u(u + t) 2 



(55) 



(56) 



where K (x) is the modified Bessel function (y = K (x) is the solution to the differential 
equation x 2 y + xy' — x 2 y = 0), so that 



dx$(x)G(xWe 



(57) 



where 



®(x) 



V2 



7T° 



2/4 K (x 2 /4) 



is a normalized distribution function ( dxQ(x) = 1). Using this to replace the inner two 
integrals of the gap equation ([521 ) we obtain 



gN(0) 



^(0) 2 J dx$(x) 



111 



2nT J V 1 2 2ttT J 



In 



D 

2^T r . 



- 7 (0) y < fa» W ^ + __j 



(58) 



where we have used the normalization of the distribution function and 7(0) 2 to extract the 
logarithm from the integral. If we set W = and T = in this expression, we obtain the 
familiar result 



1 



In 



D 



gN(0) "\2nT c0 J TX 2' 
enabling us to eliminate the coupling constant from the gap equation, writing 



1 2 ) v V 2 2ttT 



(59) 



VII. APPENDIX B 

The purpose of this section is to average the Free energy over the classical fields of the 
pseudogap. If we write the electron field in a two-component notation as 



Ck+q/2 
Ck-q/2. 
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then in the presence of the scattering potential 

H I = ^} Ac l +q /2 C ^<ll2 + A C \_ q/2 C k+q/2 ] 
k 

the corresponding electron Greens function is given by 

- (C ka {r)Cl p m = g a p(k,r) 

g a p(k,T) = T^QiK^e-^ (60) 

n 

where 

Q(k,iuj n ) = [iuj n - e k+{q/2 )r 3 - At_ - At + ]~ l 

is the inverse propagator of the electron in the pseudogap field. To obtain this expression 
we have kept scattering between states with k ~ ±k F , neglecting scattering into high energy 
states with k ~ ±(kp + nq). The free energy is then given by 



/ 



where 

F[q, A, A] = -TTrhi[-g-\k,iu; n )] (62) 

To carry out the Lorentzian integral over q, we first linearize the electron kinetic energy 
around the wavevector q = Q Q , writing 

V (g-Q ) 2 + /t 2 Trln ^" ~ e fc+(Qo/2)r 3 - v F {q - Q ) - At- - Ar + ] 

— g2 - ^, 2 Tr ]n[iw n - e k +{Q /2)r- A ~ v F q - At- - Ar + ] (63) 

This Lorentzian momentum integral has two poles in the complex plane at q — Q Q = For 
oo n > 0, the poles of the logarithm are in the upper half complex plane, so we complete the 
contour in the lower half plane, picking up the sole contribution to the pole from q = —in. 
For the opposite sign, i.e. ui n < we complete the contour in the upper half plane. The 
result of this procedure is 

F = ~ T I ^^ AA/2A2 ^H^n-e k+iQo/2)T3 -A^-Ar + ], 
u n = u n + (v F K/2)sign(u n ) (64) 
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showing that the effect of the finite correlation length is to introduce a damping term into 
the propagator. We can further simplify this result by writing the the nesting condition 

e/c+Q /2 = -efc-Q /2 as e k+Qo /2 T3 = tk-Qo/2-r?,, so that 

F = -T J ^e-^/ 2A2 Trln[^-6 fc+Qo/2 r 3 -Ar_-Ar + ], 

= ~ 2T J ^e-^Tr ln[(^ n ) 2 + {e k+Qo/2 f + AA], (65) 

To simplify the Gaussian integral over the pseudo gap amplitude, we simply make the 
change of variables 

A = V2xWe ie 

-iO 



so that 



A = V2xWe- l ° (66) 



dAdA _aa/2W 2 dxdO 



2irW 2 2tt 
so that the Free energy becomes 



F = -2T J dxe~ x Tr \n[(u n ) 2 + {t k+Qo/2 f + 2xW% (67) 
introducing a Rayleigh distribution of gap sizes. 
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